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Abstract 

Compact analytical formulae for QED radiative corrections in the processes of elastic e — p scattering are 
obtained in the case when kinematic variables are reconstructed from the recoil proton momentum measured. 
Numerical analysis is presented under kinematic conditions of current experiments at JLab. 



1 Introduction 

With the advent of new-generation electron accelerators such as CEBAF, experiments on elastic electron- 
proton scattering in the GeV-range were brought to the new level of precision. One of the novel features of 
the experiments is that in addition to the scattered electron, the momentum of recoil proton can be directly 
measured, while both the electron beam and the proton target may be polarized. The cross section depends 
only on one kinematic variable, so the measurement of final momenta gives several possibilities to reconstruct 
this kinematic variable from measurements. It results in different methods of data analysis, the best of them 
being chosen from experimental conditions such as resolution and others. Often it is used for reducing influence 
of radiative effects, which accompany any process with electrons. Ref. [El presents a review of the radiative 
corrections calculations under different approaches of kinematic variable reconstructions in experiments at 
HERA. 

In this paper we consider the case of elastic measurement 

e(fci) + P ( Pl ) — e(fc 2 ) + p(pa) (1) 

where the initial electron is polarized longitudinally and only momentum of the final proton is measured or 
is used in Q 2 = — (p% — Pi) 2 reconstruction. The corresponding method is known as reconstruction within 
hadronic variables jj], g|. There are two types of such experiments. In the first type the proton is considered 
to be polarized either longitudinally or transversely. In the second case the polarization states of recoil proton 
are measured. We calculate radiative corrections (see Fig.FJ) for these measurements and for hadronic variables 
reconstruction method. 

Measurement of final proton momentum allows to define three independent kinematic variables. A natural 
choice for them are Q 2 , azimuthal angle and so-called inelasticity u = (fci + pi — p-i) 2 — m 2 , where m is the 
electron mass. Since at the Born level u = and it > if the additional photon is emitted in the final state of 
, one can constrain (or cut) the range of u to reduce value of radiative corrections (RC). 
We calculate RC to the cross sections and asymmetries differential in Q 2 . When leptonic variables reconstruc- 
tion method is used, the integration over photonic momentum phase space cannot be performed analytically. 
The reason is that an argument of formfactors (Q 2 transferred to the proton) depends on the photon momen- 
tum. So the integration is left for numerical analysis in order to avoid additional assumptions about specific 
models for formfactors. The advantage of the reconstruction method in terms of hadronic variables is that the 
integration can be performed analytically and completely both for unpolarized and polarized parts of the cross 
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Figure 1: Feynman diagrams contributing to the Born and the radiative correction cross sections. 

section. A more detailed discussion of advantages of RC calculation for elastic scattering within leptonic and 
hadronic variables can be found in ref.|| 

There are two possibilities in considering the polarization parts of the cross section. We consider different 
polarization states. It corresponds to different coefficients of expansion of polarization vectors, while the form 
of the expansion is the same. Since these coefficients depend on inelasticity (see Section 2), which is a function 
of photon momentum, we can produce either formulae for all specific polarization cases or get only one set of 
formulae for the polarization states but leaving integration over inelasticity for numerical analysis. We choose 
the second way. We will see in Section 3, that obtained formulae are very simple even in the general case. 

We use the approach of Bardin and Shumeiko Q for this calculation. Advantages of the approach are absence 
of any approximations both in the procedure of extraction and cancellation of infrared divergent terms and in 
integration over photon phase space. However it is known that mass singularity terms still remain within these 
formulae. In our case we perform additional procedure to extract these mass singularity terms. The result of 
the expansion can be presented as 

* = log^A(Q 2 )+B(Q*)+o(^J (2) 

where independent on electron mass functions A{Q 2 ) and B(Q 2 ) appear as LO and NLO contribution to the 
cross section. The explicit formulae in the form of (|J) are presented in Section 3. 

Section 4 is devoted to numerical analysis within kinematic conditions of experiments at JLab. We analyze 
RC to unpolarized cross section and different asymmetries accessed in polarization measurements. Furthermore, 
obtained results can be applied for other reactions, where ep — > ep7 is a background process. In Section 4 a 
specific example is considered when the process is background to measurements of pion production by polarized 
real photons at JLab. 



2 The Born cross section and kinematics 

The Born cross section of the process (fD) can be written in the form 



M 2 , dQ 2 

where S = 2k\p. Kinematic limits for Q 2 are defined as 

0<Q2<^^, X s = S 2 -,m 2 M 2 (4) 

where (m, M are the electron and proton masses). Here, for generality, we keep the electron mass, however, 
below we will neglect it wherever possible. At the Born level the mass can be set to zero in the final formulae, 
while for RC cross section we will construct expansion (||). We note that in general we should keep explicit 
azimuthal angular dependence. The cross section (both Born and RC ones) does not depend on this angle, 
but kinematic cuts can. We assumed that there is only a cut on inelasticity or all other cuts can be effectively 
reduced to it. In this case integration over the azimuthal angle can be performed analytically, that results to 
(|J). Born matrix element is the contraction of leptonic and hadronic tensors: 

r2 



M o = 7Ta L % W ^ ( 5 ) 
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We use standard definitions for the leptonic tensor, 

L% = ^Tr(fc 2 + m) 7M (l + >fti)(h + m) 7 „. (6) 
The hadronic tensor can be described by four structure functions in all considered cases: 



with 



i 2 PiiPv 3 ■ wfc 4 • m fo\ 

<V = -9»u «V = -j^- = -te^xa^j-, = i.e^ A(T m3 (8) 



and (r p = Q 2 /4M 2 ) 



^ 1= 4r p M 2 G 2 n T 2 = U1 2 ^±^ F 3 = -2M*G E G M , Ta = -M 2 G M , Gm (9) 

Polarization effects are described by polarization four- vectors of the lepton (£) and proton (77) . Following || 
we expand them over the measured momenta k\ , p\ and p 2 

2 ' 5 u ^ 



N /A7 v 2m 

77 = 2a^fci +b v q + c n (px + p 2 ) (10) 

where g = p% — p 2 . In this paper we consider fourQ polarization states, when initial or final protons are polarized 
along or perpendicular (in the scattering plane) to the vector q. For these cases the following expressions of the 
coefficients are used, 

a L = 6^ Q2 + 4M2 <*= Q2 (11) 



and 



jT _ Q 2 (Q 2 +4Af 2 ) ^ bT _ (Q 2 +4M 2 )Q 2 ^ Q 2 (2S-Q 2 ) ^ 

Here Q 2 = Q 2 + u, A h = SQ 2 (S-Ql)-M 2 Qi-m 2 X M and A M = Q 2 (Q 2 + 4Af 2 ). These formulae are given for 
the case of initial particles polarizations. For the final polarizations these formulae can be kept but — > —b^. 
One can make sure that initial and final polarization vectors exactly satisfy the necessary conditions: 

Vl,t = -1) VlVt = 0, qi] T = 0, PVl,t = (13) 

where p is p\ or p 2 in dependence of initial or final polarization vector is considered. 
Calculating the contractions in (^) we obtain the Born cross section in the form 



dan 2ira 2 



dQ 2 S 2 Q 



4 E 



ObK, (14) 



where 



9° 

9 " = -^(S 2 -Q 2 S-M 2 Q 2 ) 



2Q 2 

1 
M 2 
20 2 

-^ f (Q 2 a v + (2S-Q 2 )c n ) 
9^(2S-Q 2 ){a n -2b v ) 



1 The considered model independent RC does not include box- type contributions and therefore cannot lead to the additional 
T-odd polarization and/or asymmetry. The polarization part of the cross section of the model independent RC as well as of the 
Born cross section is exactly zero for normal polarization 
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3 Radiative corrections 

The cross section of radiative process (Fig. |l|b,c) 

e(Ai)+p(pi) e'(k 2 )+ 1 (k)+p(p 2 ) : 

can be presented in a general form 



Ml 

da r = dY r 

Apk\ 



The phase space of three final particles 



1 d 3 p 2 d 3 k 2 d 3 k 



S(p + fa — k2 — k — p 2 ) 



(2tt) 5 2 P20 2k 20 2k a 

can be parameterized in terms of four invariant variables |6|, Q|: Q 2 , u = 2k 2 k, w = 2k\k and z 2 = 2p 2 k 



dT r — 



dQ 2 



8(2ir) 4 S 



du 



dw 



dz 2 

7T Z 



(15) 
(16) 

(17) 



(18) 



where R z comes from Gram determinant and coincides with a standard i? z -function from Bardin-Shumeiko 
approach ||, |(| . The formalism to calculate it is general and is developed in S . Explicit formulae for analytical 
integration are presented in papers [H [7j [| . Note that integration over the azimuthal angle is assumed to be 
performed (see discussion after Eq.(|])). 

The matrix element squared of radiated process reads 



M 2 = -p^L^W^ 

h 



The leptonic tensors of radiative process in (|T^) looks like 

1 



where 



k\ a k 2a 
kk\ kk 2 



L' 



7m 



Tr {k 2 + m)r Ma (l + 75%) (fa + m)T av , 



a 

2kki 



2kk 2 



r 

) 1 Qf 



faa k 2a 

kki kk 2 



2kk x 



2kk 2 



(19) 



(20) 



(21) 



The calculation of the cross section of the radiative process is not so straightforward because of infrared 
divergence. The separation of the infrared divergence can be performed in a standard way H : 



g?0Y = da 



IR 



da r — da. 



TR 



da. 



IR 



da,. 



The first term in r.h.s. of Eq. (|22j) has an infrared divergence and can be presented in the form: 



da 



IR 



5 IR da = -(5 s + S H )dao 



(22) 



(23) 



while the second one in (|22j) is finite for k —> 0. The quantities 63 and Sh appear after additional splitting the 
integration region over inelasticity u by the infinitesimal parameter u 



S s 



-1 



du 



du 



(27i7i)"- 4 fc 
d 3 k 



F IR 5((A h - kf - m 2 ) 



k 



F IR S((A h - kf - m 2 ) 



(24) 



where A/, = fa + p\ — p 2 and 



uw u 1 - 



w 



2 ' 



(25) 
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For this specific procedure it is convenient to keep integration over photon momentum as it is. In this case one 
delta function appears in the integrand. It is equivalent to integration over w and z 2 in (|l8|). The integration 
gives the following results 



Ss = 2[P lH + \og—){l m -l) + l + l m -l 2 m -—, 

S H = 2(/ m - 1) log - + l m ) 2 +l v + l m - l w (l m + lw - U) — -7, Li 2 ( — t^t ) , 
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Q 2 



log— l v =log-^, l w — - log(x m ), x r , 



Q 2 + Un 



(26) 



The explicit expression for the additional photon exchange contribution coming from diagram (|l|d) has also 
the factorized form: 

(27) 

= ~-2 { >" " + log ) (l m - 1) - + \hn - 2 + y (28) 



where 



day — — S dao 

7T 



Thus the sum of contributions (E3) and (E 



8 cl =S V + S S + Sjf = l m {lv - L + |) - l v ~ 1 - ~£ ~ ^2 + 2U, - Li 2 ( Im ) 



(29) 



is infrared free. 

The second term in (|22| ) is infrared free. Hence this contribution can be written using the expression for 
photon phase space (|l8|) as 



fid 



F _ a 3 dQ 2 



4 S 2 Q' 



(30) 



where only coefficients Of include integrals: 

°i =\ J du J dw J 4^=1^ ~ aFiiitfivWuv = ^ y -^-dx J dw J 
Here we use the transformation of external variable u, 



dz 



X 



Q 2 



(31) 
(32) 



to have final results in the simplest form, where x m is defined in Eq.(Pq). 

Integrals in (|3l]) can be calculated analytically completely. However, as we discussed in Introduction, one 
integration in polarization part of the cross section is retained for generality: 



9f 



9 F 



2Q 2 [—(2L -1)+L + 2l m (l w - 3) - 2l w + l 2 w - 2Li 2 (x m ) + 

1 r 7T 2 M 2 2 

J^{Q i l W ~Q 2 (S + M 2 )[L + 2l m (l w - 3) + l 2 w - 2Li 2 (z m ) + —} + -^-[1 - 2L] 

+S(Q 2 - S)[2x m (l — L) — Al w - 1] + S 2 [x 2 m L + 3 - x m - 2x 2 n - 6Z m ]} 



2Q 2 r dx 
~W J ~~x 



Q 2 a r) h 1 + (2Sh 2 + Q 2 h 3 )c v + {Q 2 Aa n + (2S - Q 2 )Ac v )h 4 



qO?F 
y 3°pol 



(33) 



91 



Q A f dx 



(Q 2 /ii + 2Sh 3 )a v + (2Sh 2 + Q 2 h 3 )b v + (2S - Q 2 )h i (Ab v - Aa v ) 
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Here the contribution of polarization part is presented in such a form where x (or u) dependence is included 
only in polarization coefficients and in functions: 



hi = {3x 2 -x + b-2l x (2x 2 + x + l))/x 2 , (34) 

hi = l-8x-2l x (l-x), 

h 3 = (5x - 5 + 2l x (x + l))/x, 

hi = (3x + i-4l x x)/(l-x). 

The special procedure of additional subtraction was applied for polarization contribution in order to extract 
mass singularity terms to l m , which is contained in l x = l m — log(a;(l — x)) and L = l m — l v + 2l w now. In 
subtracted (and added) terms arguments of polarization coefficients were taken for x — 1 (or u — 0) as for Born 
process. As a result differences A{a, b, c} n — {a, b, c} v (x) — {a, b, c}^(l) appeared in the integrands. Because of 
added terms does not contain x-dependence of polarization coefficients, they can be integrated explicitly. This 
is the origin of factorized terms with 8^ ol in expressions for 9$ and 



Spot = - I dx 



x 4 Q 2 (3Q 2 (l-x)+4xm 2 ) 



3x 3 + 3x 2 + 3x + 3 



(35) 



(Q 2 (l - x) + m 2 x) 2 

= -3(/ m + l v -l w )-l (36) 

The last contribution which has to be taken into account (see Fig. [l]e) is vacuum polarization by leptons and 
by hadrons (da vac ). Leptonic contribution comes from QED (see for example), while hadronic contribution 
can be obtained from experimental data for the process e + e~ — > hadrons. We will use a model developed in 
rcf. Q. The vacuum polarization correction has also a factorized form and we denote them S l vac and 8^ ac for 
leptonic and hadronic contributions. 

The final result for the RC cross section is obtained by adding all the considered contributions: 



4 



day + da r + da vac = £[flf + 4(S el + S l vac + (37) 



This expression is the correction to the cross section differential in Q 2 only. Sometimes it is necessary to have 
formulae for the cross section versus inelasticity also. They can be obtained by straightforward differentiation 
over x m (or u m ) of the result (p7|): 

da r da r Q 2 a 3 dQ 2 A d9? 

1^ t 38 ' 



dx du x 2 4 S 2 Q 4 ^-^ dx 

2 — 1 



where 



d0? _ n „ 2 2(l + x 2 )L + l~8x 

dx ~ Q (\-x)x 2 (M) 

d£l _ 1 ( 4(x 2 S(xS - Q 2 ) - M 2 Q 2 ) Q 4 
dx ~ M 2 1 (1 - x)x x 

(xSjxS - Q 2 ) - M 2 Q 2 )\2{\ + x 2 )L + 1 - 2x(l + a;)] ) 
(1 - x)x 2 J 
d6 R _ 2Q 2 (2(l + x 2 )L ~8x 2 + Qx-b){Q 2 a v + (2xS - Q 2 )xc v ) 
dx (1 — x)x 3 M 

d0 R _ Q 4 (2(l + x 2 )L-8x 2 + 6x-5)(2xS - Q 2 ){xb v - a v ) 
dx ~ (1 - x)x 3 M 3 

It is clear that only hard photon radiation can contribute to the functions. The terms containing x m in 8' s and 
in subtracted parts daj. R have to be completely canceled. It can be seen in polarization parts that the terms 
containing a v (l), b v (l) and c v (l) explicitly cancel out in final expressions (39). 
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4 Numerical analysis 

The obtained results can have direct application in experiments at JLab, so we perform the numerical analysis 
within kinematic conditions of JLab. First we consider the unpolarized cross section and RC to it. The next 
step is consideration of polarization effects in two cases when initial or final protons are polarized. 
Both the Born and observed cross sections can be split into unpolarized and polarized parts: 

J^bom^obs 

b.obs _i_ „b,obs I A(\\ 

dQ2 ~ " p (40) 

We consider four different polarization states described in Eqs.( pd| , |l2| ), so we have four different polarized parts 
of cross sections, which, of course, corresponds to only one unpolarized cross section. Let us define the relative 
corrections to the observable quantities in the current experiments: 



S u = 



$L,T = 



r obs 



T obs 

















U U 



f-rObS 



obs 



- 1 



5 r = 



obs 
°pT 








obs 
a pL 









T obs 
r pT 



nlis 



pL a pT 



(41) 



The first correction 5 U is the relative correction to unpolarized cross section. The Sl,t are corrections to 
polarization asymmetries measured by rotating the polarization states of initial protons. At last, the quantity 6 r 
is the correction to the measured ratio of final proton polarizations JTo| , [HJ . The correction to the unpolarized 
cross section is presented in Figure ^a. The behavior is quite typical. For the very hard inelasticity cut (u m <C 
Q 2 ) the positive contribution due to real bremsstrahlung is suppressed, so there is only negative loop correction 
contributing to cross section. Different ending values for the curves corresponds to different kinematically 
allowed regions. 

The correction to the asymmetry is given in the plot ^p. The magnitude of the correction does not exceed 
1 .5%. It goes up with increasing of Q 2 and inelasticity cut. Second effect is clear because the only unfactorized 
(different for unpolarized and polarized parts) hard bremsstrahlung can contribute to RC to asymmetries. The 
Q 2 -dependence is also understood because of contribution of a large logarithm log(Q 2 /m 2 ). Both the cross 
sections and asymmetries have non-zero leading contributions, so the correction is larger when Q 2 is going up. 
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Figure 3: Radiative correction to recoil proton polarization for (^) in the region the invariant mass of the 
unobserved state close to the pion mass. S=8GeV 2 . 
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Figure 4: The cross section and proton polarizations of the process (fl3]). Q 2 here is equivalent to —t defined 
in Q. 

In the case of calculation of correction to the transferred polarization experiment the correction is negative 



and does not exceed 1%. It is in agreement with our earlier estimation of the effect 12, 3]. 
An interesting application of our result can be found for the process 

7+p^7T°+p (42) 

recently analyzed in JLab experiment fl3|| . Usually the photons in the reaction are produced by polarized 
electron beam. As a result both photon and electrons are in the beam. Therefore when invariant mass of 
undetected e and 7 is close to the pion mass the process considered in this paper 



e+p — > e + 7 +p 



(43) 
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can be background to the process (|42|) . In the end of the section we analyze unpolarized cross section and recoil 
proton polarization due to the process (f43|). 

In Figure H^, we show the cross section vs Q 2 for several values of the unobserved invariant masses set close to 
the pion mass. Plot ||b shows polarization of protons due to the bremsstrahlung process ( [43| ) being background 
to the measurements Jl3[ . 

5 Conclusion 

In this paper we consider the radiative effects in elastic electron-proton scattering with the hadronic variables 
reconstruction method. Within this method the information on proton final momentum is used to reconstruct 
the kinematic variable Q 2 . Another kinematic variable which can be also reconstructed from the measured 
recoil momentum of proton is inelasticity. The Born cross section does not depend on the inelasticity, so this 
fact may be used in measurements to make a cut on this variables. It allows to reduce radiative corrections 
essentially, but at the cost of significant loss in statistics. In this paper we obtain explicit formulae for the 
cross section and spin asymmetries versus Q 2 and inelasticity it. We calculated RC to unpolarized cross section 
and polarization observables in kinematic conditions of experiments held at JLab. We found that correction is 
about one per cent. It increases when both Q 2 and inelasticity go up. Also we calculated a cross section and 
proton polarization from elastic ep-events that produce background in neutral pion production by polarized real 
photon at JLab. 
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